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Abstract — In this paper, a novel carrier frequency offset 
estimation approach, including preamble structure, carrier fre- 
quency offset estimation algorithm, is proposed for hexagonal 
multi-carrier transmission (HMCT) system. The closed-form 
Cramer-Rao lower bound of the proposed carrier frequency 
offset estimation scheme is given. Theoretical analyses and 
simulation results show that the proposed preamble structure and 
carrier frequency offset estimation algorithm for HMCT system 
obtains an approximation to the Cramer-Rao lower bound mean 
square error (MSE) performance over the doubly dispersive (DD) 
propagation channel. 

Index Terms — Hexagonal Multi-Carrier Transmission Sys- 
tem; Preamble Structure; Carrier Frequency Offset Estimation; 
Cramer-Rao Lower Bound; 



I. Introduction 

Orthogonal frequency division multiplexing (OFDM) sys- 
tems with guard-time interval or cyclic prefix can prevent 
inter-symbol interference (ISI). OFDM has overlapping spectra 
and rectangular impulse responses. Consequently, each OFDM 
sub-channel exhibits a sinc-shape frequency response. There- 
fore, the time variations of the channel during one OFDM sym- 
bol duration destroy the orthogonality of different subcarriers, 
and result in power leakage among subcarriers, known as inter- 
carrier interference (ICI), which causes degradation in system 
performance. In order to overcome the above drawbacks of 
OFDM system, several pulse-shaping OFDM systems were 
proposed HMD. 

It is shown that signal transmission through a rectangular 
lattice is suboptimal for doubly dispersive (DD) channel (9)- 
[ 1 3 1 . By using results from sphere covering theory, the au- 
thors have demonstrated that hexagonal multi-carrier transmis- 
sion (HMCT) system obtains lower energy perturbation than 
OFDM system, hence providing better performance against 
ISI/ICI 0-03]]. There are mainly two types of HMCT 
systems, (1): HMCT system with orthogonal prototype pulse 
(Hl-QO), named as lattice OFDM (LOFDM) system, which 
confines the transmission pulses to a set of orthogonal ones. 
(2): HMCT system with nonorthogonal prototype pulse [ 1 1 1- 
[13 1, named as hexagonal multicarrier transmission (HMT) 
system, which abandons the orthogonality condition of the 
modulated pulses and obtains the high robustness performance 
on combating the dispersion of DD propagation channel in 
both time and frequency domain. 



To optimally combat the impact of the DD propagation 
channels, the lattice parameters in the time-frequency (T- 
F) plane and the pulse shape of modulation waveform of 
HMCT system are jointly optimized to adapt to the channel 
scattering function. It is shown in |14|-[17| that the HMCT 
system outperforms OFDM system from the robustness against 
channel dispersion point of view. 

The basic mathematical operation of the received signal 
performed by the demodulator is a projection onto an identi- 
cally structured function set generated by the prototype pulse 
function [5], i.e. an optimal match filter. In HMCT system, 
there is no cyclic prefix and data symbols of HMCT signal 
are transmitted on hexagonal lattice points in the T-F plane. 
Traditional timing and frequency synchronization schemes can 
not be applied to HMCT system directly. In our previous work 
lfT4l - lfT7l . signal to interference and noise ratio (SINR) of 
HMCT system over DD channel is analyzed and the Max- 
SINR based timing synchronization scheme is proposed. 

In this paper, a novel preamble structure and two carrier 
frequency offset estimation algorithms, named as phase dif- 
ferential algorithm and least squares algorithm, for HMCT 
system are proposed. After detailed derivation, the close form 
performance lower bound of the proposed carrier frequency 
offset estimation algorithm is given. Theoretical analyses and 
simulation results show that the proposed carrier frequency 
estimation approach for HMCT system obtains an approxima- 
tion to the lower bound mean square error (MSE) performance 
over DD channel. 

II. Hexagonal Multi-Carrier Transmission 
System 

In HMCT systems, the transmitted baseband signal can be 
expressed as iTTTl — lfT3l 
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where T and F are the lattice parameters, which can be 
viewed as the symbol period and the subcarrier separation, 
respectively; c TO ,„ denotes the transmitted data symbol, which 
is assumed to be taken from a specific signal constellation 
and independent and identically distributed (i.i.d.) with zero 



mean and average power a 2 ; m G M and n G N are 
the position indices in the T-F plane; M. and M denote 
the sets from which m, n can be taken, with cardinalities 
M and N, respectively. ip m ,2n(t) = - mT)e ,2l " Fi and 
ipm,2n+i(t) = ip(t-mT-^)e^ 7Ti - nF+ ^ t are the transmitted 
pulses generated by the prototype pulse ijj(t). The ambiguity 
function of the prototype pulse is defined as 
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-j2™t dt 
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The baseband DD channel can be modeled as a random linear 
operator H [18] 



U[x(t)} 



fd 



H(t, v)x{t - T)e j27TVt drdv (3) 



-fd 



where T max and fd are the maximum multipath delay spread 
and the maximum Doppler frequency, respectively [19|. 
H(t,v) is called the delay-Doppler spread function, which 
is the Fourier transform of the time-varying impulse response 
of the channel h(t, t) with respect to t Ifflfl . 

In wide-sense stationary uncorrelated scattering (WSSUS) 
assumption, the DD channel is characterized by the second- 
order statistics 



E[H{t, v)H*(tu Vi )] = S h (t, v)S{t - n )S(t 
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(4) 



where E[-] denotes the expectation and Sh(t,v) is called 
the scattering function, which characterizes the statis- 
tics of the WSSUS channel. Without loss of generality, 
H(t,v) is assumed to have zero mean and unit variance, 
i.e.,J^' m J^ d j^SH(T,v)dTdv = 1. The received baseband 
signal can be expressed as 



r(t) =U[x(t)] +w{t) 



(5) 



where w(t) is the AWGN with variance crj. 

The basic mathematical operation of the received signal 
performed by the demodulator is a projection onto an identi- 
cally structured function set generated by the prototype pulse 
function, i.e. an optimal match filter 1171 . To recover the trans- 
mitted data symbol c m ,2n, the match filter receiver projects 
the received signal r(t) on prototype pulse function i/jm,2n(t), 

i.e., c m>2n = (r(t),ip m ,2 n (t)) = /!^ o r(t)^ i 2n(*) d *> where 
(•) denotes the inner product and (•)* denotes the complex 
conjugate. Hence, frequency synchronization plays a critical 
role in ensuring reliable demodulation. Firstly, we will analyze 
the effects of carrier frequency offset on the recovered data 
symbol c m , 2n . 

III. Effects of Carrier Frequency Offset on HMCT 
System 

Under the assumption that there is a carrier frequency offset 
A/ between the received signal r(t) and the transmitted signal 
x(t). After ignoring the impact of additive noise, the received 
signal r(t) can be written as 



r(t) 



= e j27rA/t 



fd 



H(t, v)x(t - T)e j27TVt dTdv (6) 



and the recovered data symbol c m: 2n can be expressed as 

C m ,2n = (r(t),ljj m ,2n(t)) 
N/2-1 



Cm ' , 2 «' ,£1 mfn;m' ,2ri 
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(7) 
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where S^ n . m , 2n , in (7) can be expressed as (8) at the top of 
the next page. 

Let m = m' and n — n\ S^n-m 2n can be expressed 



' — 'm,2n;m,2n 



and H m „. m , 2n '+i m (7) can ^ e expressed as (10) at the top 
of the next page. 

Let m = m! and n — n\ S m 2 n -m 2n+i can ^e expressed 

as 
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■ e^ v( - m+ ^ T drdv 
Hence, the data symbol c m ^n in (7) can be rewritten as 



_ -A/ \ " _A/ 

C m ,2n — c m,2n^ m; 2n;m,2n + 2-^1 Cm , 2n ' "m,2n;m,2n 
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The first term in equation (12) represents the desired symbol 
and the last three terms denote ISI/ICI. Concerning the useful 
portion, the transmitted symbols c m> 2n are attenuated by 
j4ff(7inax> fd, A/) which is caused by the carrier frequency 
offset A/ and doubly dispersive channel. Meanwhile, the 
transmitted symbols rotated by a phasor —j2irmTAf. 

IV. The Proposed Preamble Structure 

The proposed preamble is composed of two training se- 
quences Pi and P-2 in the frequency domain, as depicted in 
Fig. 1. P, = [Pi(0),Pi(l),--- , Pi(N P — 1)], i G {1,2} and 
Np < N/2 denotes the length of training sequence. 

In this paper, Pi are selected as PN sequence and Pi = 
P2. Assume that the index of the training sequence is I G 



TtA/ II I J2-wAft 
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Fig. 1. The proposed preamble structure. 

[0, 2Np — 1], the frequency domain training sequence can be 
expressed as 

(Pi([l/2\), l€[0,2N P -l] 
° U = \ , else (13) 

where |_-J denotes the floor function. 

V. The Proposed Carrier Frequency Offset 
Estimation Algorithm 

We assume that the received signal r(t) experiences a 
carrier frequency offset Af and the timing offset is completely 
compensated. The demodulated training sequences can be 
expressed as 

N/2-1 

Dl,2l=Dl,2l^t t 2l;l,2l + E D UI ,S UI;1,2I' 
l'=0,l'jil 

N/2-1 
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k',2l'+l^f\,2l+l;k',2l' + l 
k'jtl l'=0 
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Wk 21 ar, d Wk,2i+i denote the AWGN in frequency domain. 
Thanks to the central limit theorem ll22l . the last four terms in 
(14) and (15) can be denoted as 4>\.2i and 4>i^i+i, respectively. 
We can rewrite equation (14) and (15) as 



Dl,2l = £>l,2/Sf^ ;1 2 , + 01,2/ 

= n..2,- - I Af 

■ e-^ 2lF ' T e^ vT 'drdv 



h(r,v)A* i ,(T,v + Af) (16) 



>l,2l 



and 



Diai+l — D 1,21+1 ^l,2l+l;l,2l+l + $1,21+1 

= ^, 2 me- j27r(1+ ' )T ' A/ / / Hr,v) 

. ^1+1/2)1^ dTdv + (h 



(17) 



We assume that Di 21D* 2 ; — u s and define ^y m as (18) 
at the top of next page. $0 = irT'Af denotes the phase 
rotation caused by frequency offset and 1J1 is a constant. 
Let e = AfT s N/2 denotes the normalized frequency offset, 
and T s denotes the sample period. W is the additive noise 
with zero mean and average power cr 2 ^. There are two 
fine frequency offset estimation schemes: Phase differential 
estimation scheme and least squares estimation scheme. 



N/2-1 . . * 
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A. Phase Differential Carrier Frequency Offset Estimation 
Algorithm 

The phase differential carrier frequency offset estimation 
algorithm can be expressed as 



N 
2ttM 



arg(7o) 



; N/2 - 1 \ 
N /2-l E ( ar S(7m) -arg(7m-i))J 



(19) 



The phase differential estimation scheme is valid for e e 

[-N/2M,N/2M] = [-p/2,p/2], and p denotes system 
spectral efficiency. 

B. Least Squares Carrier Frequency Offset Estimation Algo- 
rithm 

We can model the phase of -f m as 

arg( 7m ) ~ i? + (m + + Cm (20) 

and 

arg(7) ~ CS+ C (21) 

where 7 = [70,71, • • • ,7w/ 2 -i] ^ and arg(7) = 
[arg( 7o ),arg( 7 i),--- ,arg( 7A r /2 _i)] T , # = [i? ,#i] and 
C = [Co, Ci, • • • , Cjv/2-i]- C in (21) can be expressed as 

1 



c = 



1 



1 



(22) 



.1 ^72-1. 

The least squares carrier frequency offset estimation algorithm 
can be expressed as 



= (C T C)- 1 C T arg( 7 ) 
and (C T C) _1 in (23) can be obtained by 



(23) 



(C^C)- 1 



NAm 2 



m z — m 
—m 1 



(24) 



where m 



= l/N^ll 1 * 2 and 



Am 2 = 2/Nj2k=l 1 (k ~ m ) 2 - Hence, the estimation of $ 
can be rewritten as 

1 



0= [<Mif 



Am 2 



7m 2 m(7m) 
Tto m • 7 



(25) 



where 7 = 2/TV^^ 2 , 1 arg(7^) and 7m = 

2/NJ2k=o 1 aT S(lkk). The phase rotation #0 caused 
by frequency offset and the constant #1 can be obtained by 



01 



7m — m • 7 
ZW 2 



and 



<9 = arg(7o) - i?i = 



7m 2 — 771(777!,) 
Am 2 



(26) 



(27) 



The normalized frequency offset e can be expressed as 

= iV(arg(7o) -tfi) (2g) 
2?rM 

The least squares frequency offset estimation scheme is valid 

for e e [-N/2M, N/2M] = [-p/2,p/2]. 

C. Simplified Carrier Frequency Offset Estimation Scheme 

Let * = hi>o(0),Vo(l),-- - ,Vo(iv ~ !)] denotes the 
prototype pulse with length L$ in discrete HMCT system. 
In order to recover the transmitted frequency domain training 
sequence D\ t i, we need N p inner product operations with L^, 
points, which introduce high computational complexity. 

Theorem 1: Let r = [r(0), r(l), • • • , r(L^ — 1)] denote the 
received signal vector with length L^, and <& = {*f? n }> n = 
0, 1, • ■ • , N/2 — 1, denotes the objective projection subspace. 
*n = [ipn{0),ipn(l),--- ,i>n{L^ - 1)] and satisfies 



ip n (m) = ip {m)e" 



JV/2 



(29) 



The time frequency subspace projection of the received 
signal r on the objective projection subspace <& is equivalent 
toFFT(Y / ?J ~ 1 Mk + lN/2)), k ^ 0,1, - ■ ■ , N/2-1. = 
r^/sl anc * PI denotes the ceiling function. r^rQfJ and 
r$ = [t>(0), t\/,(1), • • ■ , r^,(L^ — l)]. denotes the Hadamard 
product and FFT(-) denotes the Fast Fourier Transform. 

Proof: The output symbol Y(i) of the time frequency 
subspace projector after projecting the received signal r to the 
ith subspace \&i can be expressed as 

Lw, — 1 



K(») = (r,* < >= ]T r{n)i)*(n) 



n=0 



L*b — 1 



(30) 



r(n)^(n)e ~/ 2 
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After performing the IFFT on the signal Y (i), the transformed 
signal can be expressed as 



Differentiating the log likelihood function lnP( 7 o;e) with 
respect to e, we have 



j27Tki 
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(31) 



i=0 n=0 ;=0 



£jV\ j2^(fc-n)i 
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E r ^( fc +-o-) 



where ["•] denotes the ceiling function. r^,(n) = r(n)f/'o( n )> 
n = 0, 1, • • • , L$ — 1, which can be denoted as = r0 *q. 

■ 

We can conclude from Theorem 1 that the transmitted 
frequency domain training sequence Dij can be obtained by 
the following two steps. The first step is to calculate the 
Hadamard product of the received signal r and the prototype 
pulse *o, that is r$. The second step is to superimpose 
vector with period N /2 and perform 7V/2-point FFT on 
the superimposed vector y = [y(0), 2/(1), y(N/2)]. Hence, 
two Hadamard product operations, two superimpose operators 
with period N/2 and two 7V/2-point FFT are needed to 
recover the transmitted frequency domain training sequence 
D\ t i, which is a low complexity approach compared to the 
traditional N p projectors approach. 

D. Cramer-Rao Lower Bound 

We can conclude from (19) and (28) that the proposed 
phase differential and least squares carrier frequency offset 
estimators are functions of 70. 70 in AWGN channel can be 
expressed as 

N/2-1 t 

7o= E Ue- jvT ' Af + ^,2i + ia s e^ T '^ 
1=0 { 



(32) 



N 



a 2 e -o2nMe/N + Wq 



The conditional probability density function P(7o;e) can be 
expressed as 

_. . 1 f 1 / N „ -,2.M. \ 

P (7o; e) = 7—2 u /2 ex P1 " Z2— (To - y^ s e » J 

7TOW ' ^ OW- V 2 / 



,7o-yf- s e « 



■ { - j7rMa 2 e^J + (70 - — a^e" 



(34) 



and 



a 2 lnP( 7 o;e) 1 /2tt 2 M 2 ct 2 22^ 
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Hence, the Cramer-Rao lower bound of the proposed carrier 
frequency offset can be expressed as 



E{|e-e| 2 } > -E 



a 2 lnP( 70 ;e) 



d 2 e 



N 



(36) 



(33) 



2tt 2 M 2 SNR 
VI. Simulation Results 

In this section, we test the proposed synchronization ap- 
proach for MCM system with hexagonal T-F lattice via com- 
puter simulations based on the discrete signal model. In the 
following simulations, the number of subcarriers for HMCT 
system is chosen as iV=40, and the length of prototype pulse 
L,/,=600. The center carrier frequency is / c =5GHz and the 
sampling interval is set to T„=10 _6 s. The system parameters 
of HMCT system are P=25kHz, T=l x l(T 4 s and signaling 
efficiency is set to p=0.8. WSSUS channel is chosen as DD 
channel with exponential power delay profile and U-shape 
Doppler spectrum. 

The MSE performance of the proposed carrier frequency 
estimation algorithm over AWGN channel is given in Fig. 2. 
We can see from Fig. 2 that the proposed least squares carrier 
frequency estimation algorithm outperforms phase differential 
scheme at low SNR, and both the proposed schemes can 
obtain an approximation to the Cramer-Rao lower bound MSE 
performance. 

The MSE performance of the proposed carrier frequency 
estimation algorithm over DD channel is given in Fig. 3. 
As shown in Fig. 3 that the proposed least squares carrier 
frequency estimation algorithm outperforms phase differen- 
tial scheme at low SNR, but the phase differential scheme 
outperforms least squares algorithm at high SNR. The DD 
propagation channel introduces energy perturbation among 
the transmitted symbols, hence there is a gap between the 
MSE performance of the proposed two carrier frequency 
offset estimation algorithms and that of the Cramer-Rao lower 
bound. 
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Fig. 2. MSE performance of the proposed carrier frequency estimation 
algorithm over AWGN channel. 
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Fig. 3. MSE performance of the proposed carrier frequency estimation 
algorithm over DD channel. 



VII. Conclusion 

In this paper, the effects of carrier frequency offset on 
HMCT system is analyzed. Meanwhile, a novel carrier fre- 
quency offset estimation approach, including preamble struc- 
ture, carrier frequency offset estimation algorithm, is proposed 
for HMCT system. After detailed derivation, the close form 
performance lower bound of the proposed carrier frequency 
offset estimation alglrithm is given. Theoretical analyses and 
simulation results show that the proposed carrier frequency 
offset estimation approach for HMCT system obtains an 
approximation to the Cramer-Rao lower bound MSE perfor- 
mance over the DD channel with exponential power delay 
profile and U-shape Doppler spectrum. 
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